Kinetic properties of particle-in-cell simulations compromised by Monte Carlo collisions

he particle-in-cell method with Monte Carlo collisions is frequently used when a detailed kinetic simulation of a weakly collisional plasma is required. In such cases, one usually desires, inter alia, an accurate calculation of the particle distribution functions in velocity space. However, velocity space diffusion affects most, perhaps all, kinetic simulations to some degree, leading to numerical thermalization (i.e., relaxation of the velocity distribution toward a Maxwellian), and consequently distortion of the true velocity distribution functions, among other undesirable effects. The rate of such thermalization can be considered a figure of merit for kinetic simulations. This article shows that, contrary to previous assumption, the addition of Monte Carlo collisions to a one-dimensional particle-in-cell simulation seriously degrades certain properties of the simulation. In particular, the thermalization time can be reduced by as much as three orders of magnitude. This effect makes obtaining strictly converged simulation results difficult in many cases of practical interest

Nonlinear resonance absorption in laser-cluster interaction

Rare gas or metal clusters are known to absorb laser energy very efficiently. Upon cluster expansion the Mie plasma frequency may become equal to the laser frequency. This linear resonance has been well studied both experimentally and theoretically employing pump probe schemes. In this work we focus on the few-cycle regime or the early stage of the cluster dynamics where linear resonance is not met but nevertheless efficient absorption of laser energy persists. By retrieving time-dependent oscillator frequencies from particle-in-cell simulation results, we show that nonlinear resonance is the dominant mechanism behind outer ionization and energy absorption in near infrared laser-driven clusters.Comment: 4 pages, 4 figures, REVTeX, minor modifications according to referee comments, accepted for publication in Phys. Rev. Let

Limiting stable currents in bounded electron and ion streams

The classical static analysis of the infinite planar diode has been extended to include the effects of finite transverse beam size. Simple expressions have been found for the increase in maximum stable current density over that of an infinite stream for finite cylindrical and strip streams flowing between plates of infinite diodes. The results are also given in terms of stream perveance. The effect of a nonuniform distribution of current across the stream is shown to be relatively small. Experimental values of maximum stable current agree with those obtained from the analysis. A further extension of the static analysis has been made to include the effects of additional conducting plane boundaries parallel to the stream motion. For length-to-width ratios L/D less than 0.25 the tube is adequately described by the results for the infinite planar diode and for L/D greater than 4, the infinitely-long drift tube theory suffices. At intermediate values of L/D, the maximum amount of current that can be stably passed through the tube is greater than that predicted by either asymptotic theory

Collisionless heating in capacitive discharges enhanced by dual-frequency excitation

We discuss collisionless electron heating in capacitive discharges excited by a combination of two disparate frequencies. By developing an analytical model, we find, contrary to expectation, that the net heating in this case is much larger than the sum of the effects occurring when the two frequencies act separately. This prediction is substantiated by kinetic simulations, which are also in excellent general quantitative agreement with the model for discharge parameters that are typical of recent experiments

Effects of high energy photon emissions in laser generated ultra-relativistic plasmas: real-time synchrotron simulations

We model the emission of high energy photons due to relativistic charged particle motion in intense laser-plasma interactions. This is done within a particle-in-cell code, for which high frequency radiation normally cannot be resolved due to finite time steps and grid size. A simple expression for the synchrotron radiation spectra is used together with a Monte-Carlo method for the emittance. We extend previous work by allowing for arbitrary fields, considering the particles to be in instantaneous circular motion due to an effective magnetic field. Furthermore we implement noise reduction techniques and present validity estimates of the method. Finally, we perform a rigorous comparison to the mechanism of radiation reaction, and find the emitted energy to be in excellent agreement with the losses calculated using radiation reaction

Inductive and Electrostatic Acceleration in Relativistic Jet-Plasma Interactions

We report on the observation of rapid particle acceleration in numerical simulations of relativistic jet-plasma interactions and discuss the underlying mechanisms. The dynamics of a charge-neutral, narrow, electron-positron jet propagating through an unmagnetized electron-ion plasma was investigated using a three-dimensional, electromagnetic, particle-in-cell computer code. The interaction excited magnetic filamentation as well as electrostatic plasma instabilities. In some cases, the longitudinal electric fields generated inductively and electrostatically reached the cold plasma wave-breaking limit, and the longitudinal momentum of about half the positrons increased by 50% with a maximum gain exceeding a factor of 2 during the simulation period. Particle acceleration via these mechanisms occurred when the criteria for Weibel instability were satisfied.Comment: Revised for Phys. Rev. Lett. Please see publised version for best graphic

Power loss in open cavity diodes and a modified Child Langmuir Law

Diodes used in most high power devices are inherently open. It is shown that under such circumstances, there is a loss of electromagnetic radiation leading to a lower critical current as compared to closed diodes. The power loss can be incorporated in the standard Child-Langmuir framework by introducing an effective potential. The modified Child-Langmuir law can be used to predict the maximum power loss for a given plate separation and potential difference as well as the maximum transmitted current for this power loss. The effectiveness of the theory is tested numerically.Comment: revtex4, 11 figure

Ultraviolet avalanche in anisotropic non-Abelian plasmas

We present solutions of coupled particle-field evolution in classical U(1) and SU(2) gauge theories in real time on three-dimensional lattices. For strongly anisotropic particle momentum distributions, we find qualitatively different behavior for the two theories when the field strength is high enough that non-Abelian self-interactions matter for SU(2). It appears that the energy drained by a Weibel-like plasma instability from the particles does not build up exponentially in transverse magnetic fields but instead returns, isotropically, to the hard scale via a rapid avalanche into the ultraviolet.Comment: 22 pages, 10 figures; v3: small textual changes; updated to correspond with version to appear in publicatio

X-band microwave generation caused by plasma-sheath instability

It is well known that oscillations at the electron plasma frequency may appear due to instability of the plasma sheath near a positively biased electrode immersed in plasma. This instability is caused by transit-time effects when electrons, collected by this electrode, pass through the sheath. Such oscillations appear as low-power short spikes due to additional ionization of a neutral gas in the electrode vicinity. Herein we present first results obtained when the additional ionization was eliminated. We succeeded to prolong the oscillations during the whole time a positive bias was applied to the electrode. These oscillations could be obtained at much higher frequency than previously reported (tens of GHz compared to few hundreds of MHz) and power of tens of mW. These results in combination with presented theoretical estimations may be useful, e.g., for plasma diagnostics.Comment: 12 pages, 7 figure

Nonlinear electrostatic oscillations in a cold magnetized electron-positron plasma

We study the spatio-temporal evolution of the nonlinear electrostatic oscillations in a cold magnetized electron-positron (e-p) plasma using both analytics and simulations. Using a perturbative method we demonstrate that the nonlinear solutions change significantly when a pure electrostatic mode is excited at the linear level instead of a mixed upper-hybrid and zero-frequency mode that is considered in a recent study. The pure electrostatic oscillations undergo phase mixing nonlinearly. However, the presence of the magnetic field significantly delays the phase-mixing compared to that observed in the corresponding unmagnetized plasma. Using 1D PIC simulations we then analyze the damping of the primary modes of the pure oscillations in detail and infer the dependence of the phase-mixing time on the magnetic field and the amplitude of the oscillations. The results are remarkably different from those found for the mixed upper-hybrid mode mentioned above. Exploiting the symmetry of the e-p plasma we then explain a generalized symmetry of our non-linear solutions. The symmetry allows us to construct a unique nonlinear solution up to the second order which does not show any signature of phase mixing but results in a nonlinear wave traveling at upper-hybrid frequency. Our investigations have relevance for laboratory/astrophysical e-p plasmas